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We propose an experiment in which long wavelength discrete axial quasiparticle modes can be 
imprinted in a 3D cigar-shaped Bose-Einstein condensate by using two-photon Bragg scattering 
experiments, similar to the experiment at the Weizmann Institute [J. Steinhauer et al, Phys. Rev. 
Lett. 90, 060404 (2003)] where short wavelength axial phonons with different number of radial 
modes have been observed. We provide values of the momentum, energy and time duration of the 
two-photon Bragg pulse and also the two-body interaction strength which are needed in the Bragg 
scattering experiments in order to observe the long wavelength discrete axial modes. These discrete 
axial modes can be observed when the system is dilute and the time duration of the Bragg pulse is 
long enough. 

PACS numbers: PACS numbers: 03.75.Kk, 32.80.Lg, 67.40.Db 



Bose-Einstein condensates (BEC) [l| of alkali-atoms 
represent a wonderful testing ground of theories of weakly 
interacting bosons. The atoms are spatially confined by 
a harmonic trap and inhomogeneity makes these systems 
even more interesting. By changing the trapping fre- 
quency one can produce quasi-lD, quasi-2D as well as 
cigar-shaped but 3D Bose systems. Due to the reduc- 
tion in the dimension and different shape, the effect of 
density and phase fluctuations become prominent which 
produces many interesting features. Bragg spectroscopy 
[3, El of a trapped BEC has become an important tool to 
reveal many bulk properties such as dynamic structure 
factor, verification of the Bogoliubov excitation spectrum 
0, 0| , observation of the Bogoliubov quasiparticle ampli- 
tudes [||, momentum distribution and correlation func- 
tions of a phase fluctuating quasi- ID Bose systems 
Bragg spectroscopy has also been used for tuning and 
measuring the scattering length of atoms by using the 
optical Feshbach resonance [9(- 

In the Bragg scattering experiments, the condensate 
is excited by using two Bragg pulses with approximately 
parallel polarization, separated by an angle 9. The pulses 
have a frequency difference u> determined by two acousto- 
optic modulators. If a photon is absorbed from the 
higher- frequency (wh) beam and emitted into other (u>i), 
an excitation is produced with energy fiuj = h(u>h — loi) 
and momentum Kk, where k = |k| = 2k p sin(9/2), and 
k p is the photon wave number. The wave-vector k is ad- 
justed to be along the z-axis. Both the values of k and lj 
can be tuned by changing the angle between two beams 
and varying their frequency difference. For small values 
of fc, the system is excited in the phonon regime and the 
response is detected by measuring the net momentum 
imparted to the system. Note that, all the momentum 
response experiments 0- [J are limited to at most a 
quarter of radial trapping period To (To = 2ir/uJo) and 
the energy transfer u) is chosen to be the order of kHz so 
that u> > too (where loq is the radial trapping frequency) 



which excites the different radial modes. In those exper- 
iments 0, 0] , the local density approximation has been 
used to the actual inhomogeneous condensates, by using 
the Bogoliubov theory of uniform gases. 

The axial excitations of a cigar shaped BEC can be 
divided into two regimes: short wavelength excitations 
whose wavelength is much smaller than the axial size and 
long wavelength excitations whose wavelength is equal 
or larger than the axial size of the system. The short 
wavelength axial phonons with different number of ra- 
dial modes of a cigar-shaped condensates which give rise 
to the multibranch spectrum has been resolved in a 
Bragg spectroscopy with a long duration (is > T ) of the 
Bragg pulses 0, 03 ■ In this experiment the condensate 
is excited by the Bragg pulses and allow it for free ex- 
pansion. During the expansion of the condensate, these 
short wavelength excitations decouple from the conden- 
sate itself, forming a separate cloud of atoms. The total 
momentum is measured by counting the atoms which is 
decoupled from the condensate. 

We consider a 3D Bose gas confined in a cigar-shaped 
harmonic trap. As described in Ref. 0|, the conden- 
sate consists of N = 10 5 atoms of 87 Rb. The radial 
and axial trapping frequencies are loq = 2n x 220 Hz 
and w z = 2ir x 25 Hz (a z = ^Jh/mu) z = 2.155/ito), re- 
spectively. Recently, the low-energy axial modes of a 3D 
BEC has been studied theoretically and experimen- 
tally Q by using a matter-wave interference technique. 
In those works [la. flij , they have shown that the long 
wavelength phonon do not form a separate cloud dur- 
ing the expansion. At long expansion times and for a 
given k, the excitations can separate out from the con- 
densate only if kag > k c a — n(u z /U( ) )(r]/2) 1 i 2 , while for 
k < k c , they remain within the condensate at all times 
|l3|. Here, r\ = [i/Tiujq and /1 is the chemical potential. 
Therefore, in order to observe the long wavelength dis- 
crete axial modes by using the Bragg spectroscopy as in 
Ref. , we have to choose the parameters such that the 
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low energy excitations can separate out from the con- 
densate for low value of k. This is an alternative way to 
realize the low-energy discrete axial modes in contrast to 
the matter- wave interference technique [l4| . 

In order to get the separate cloud of low-energy axial 
modes, the chemical potential r\ must be small. One can 
reduce the two-body scattering length a by using the Fes- 
hbach resonance [2J so that the axial size decreases and 
the long wavelength excitations can separate out easily 
from the condensate. If we take a — 13.63 x 10 _11 m, 
instead of the bare scattering length for Rubidium and 
other parameters are remain fixed as it is in Ref. ||, 
then fi — 2hu>Q. Of course, one can also get the chemical 
potential fi — 2%ujq by reducing the number of atoms in 
the condensate while keeping the bare value of the two- 
body scattering length. For this choice of a, k c ao = 0.35 
and we will show that the long wavelength excitations 
can be separated out of the condensate. Note that for 
this choice of scattering length, the system is still be- 
ing 3D. The radial and axial sizes of the condensate are 
Rq ~ 1.46/im and Zq ~ 12.78/xto, respectively, so that it 
becomes a cigar-shaped condensate with the deformation 
parameter A = luq/u> z ~ 9. 

Note that in the experiment 0, the energy transfer u> 
is chosen to be the order of kHz (ui > qjq) which excites 
the short wavelength axial phonons with different radial 
modes. Since the wave-lengths of these excitations are 
very small compared to the axial size and large compared 
to the radial size, then the axial quasiparticle modes can 
be described by the plane-wave states and the radial ex- 
citations can be described by the discrete modes which 
are the manifestation of the inhomogeneity along the ra- 
dial direction. Therefore, the short wavelength phonons 
propagate along the z-axis with different number of ra- 
dial modes which have been observed in Ref. If oj 
is the order of uj z but less than the radial trapping fre- 
quency ujq i.e. lo z < lo « wo, and the wave- vector of 
the Bragg pulse is comparable to the inverse of the axial 
size then it would excite only the different number of long 
wavelength axial modes, instead of short wavelength axial 
phonons with different numbers of radial modes. There- 
fore, the wave-length of these excitations are comparable 
to the axial size and the effect of inhomogeneity along 
the axial direction has to be included in the studies. The 
inhomogeneity along the axial direction would be mani- 
fested by the presence of the discrete axial modes. The 
discrete long wavelength axial modes due to the finite 
size of the axial direction can be observed by measuring 
the dynamic structure factor, S(k, ui), which is related to 
the momentum transferred P z (t) due to the low-energy 
two-photon Bragg scattering. 

Within the Thomas-Fermi approximation ( i. e. ne- 
glecting the quantum pressure term), the hydrodynamic 
description for the axial density fluctuations Sn(z) of a 
3D cigar shaped Bose system has been studied in |l5j . In 
the dimensionless form, the equation for the 5n{z) reads 



as 

[(1 - 5 2 )V? - 4zV 2 ~ + ACj 2 ]5n{z ) = 0, (1) 

where z = z/Zq and Co — u)/ui z . The eigenvalues of the 
above equation are {ujj/uj z ) 2 = j(J + 3)/4 and the corre- 
sponding normalized eigen functions are 

where P^ X,1 \z) is the Jacobi polynomial. Eq. (JJJ, in- 
cluding the above eigenvalues and the eigen vectors are 
valid only when the system parameter satisfies the con- 
dition: fj, » hujQ » huj z . In our method, the chemical 
potential is not too large compared to the first radial ex- 
citation, but it is quite large compared to the first axial 
excitation. Therefore, the Thomas-Fermi approximation 
is a reasonable approximation, since we are dealing with 
the axial modes only and there is no coupling with the 
radial modes. The quantum numbers j = 1 and j = 2 
correspond to the axial center-of-mass and the breathing 
modes with the frequencies u>\ = u> z and uii — \fh/2u> z , 
respectively. Only these two frequencies have been mea- 
sured with high accuracy by using the time-dependent 
modulation of the trapping potential \\(\ . It is very diffi- 
cult to measure the frequencies of the other axial modes 
of large quantum numbers by using the time-dependent 
modulation of the trapping potential and therefore there 
are no measurements of the frequencies of these higher 
modes. It is useful to verify the frequencies of the other 
axial modes (J > 2) to make sure that the energy eigen 
values derived from the hydrodynamic approximation are 
correct. These modes can be easily verified by using the 
Bragg spectroscopy which is also our main concern of this 
work. 

The dynamic structure factor is obtained from the 
Fourier transform of the time-dependent density-density 
correlation functions, 

S(k,u) = J dtj dze l{ut - kz) < 5h(z,t)Sn(0,0) >, (3) 

where the density fluctuation operator is 



6h(z,t)=J2i (4) 

Here, g = AiraTi 2 /m is the interaction strength deter- 
mined by the s-wave scattering length a and a.j is the 
quasi-particle operator of the j-th mode. It is the den- 
sity fluctuation spectrum that can be measured in the 
two-photon Bragg spectroscopy. At T = 0, the dynamic 
structure factor can be written as 

S(k,u) = Yl A Mi( k )\ 2s ("-^)> ( 5 ) 
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where A; 



Pj 14) (5). Here, k 



cast for this system as, 



P!- ' (z) is the spatial Fourier transform of 

kZo is the dimensionless wave 
vector. We rewrite the dynamic structure factor as 
S(k,u) = Ej - Uj), where Sj(k) = A,-|^-(fc)| 2 

is the weight factor which determines the weight of the 
Bragg-scattering cross-section in S(k,uj) of the corre- 
sponding axial modes of energy TwJj. In Fig.l we show 
Sj(k) as a function of the dimensionless wave vector 
k = kZo for the excitations j — 1, 2, and 3. Fig.l shows 
how many axial modes significantly contribute to S(k, ui) 
for a given momentum fik. It is clear from Fig.l that 
the strongest weights for these collective modes appear 
for kZo > 2. As an example, for our choice of parame- 




FIG. 1: (Color online) Plots of the weight factor Sj(k) vs the 
dimensionless wave-vector k for j — 1, 2, and 3. 

ters which gives Z$ ~ 5.932a z , this means that the mo- 
mentum transfer Tik in a Bragg scattering experiments 
should be k > 0.337OJ 1 for N — 10 5 in order to pick up 
the strong spectral weight from the low-energy collective 
modes. However, in order to get the separate cloud of 
the low-energy excited atoms from the condensate, the 
wave- vector of the Bragg pulse must be k > 0.35<2q 1 
which implies that the wave vector of the Bragg pulse 
must satisfy the condition: kZo > 6. 

We compare the dynamic structure factor given in Eq. 
(jSJ) with the dynamic structure factor calculated within 
the local-density approximation (LDA) . In the LDA, the 
dynamic structure factor is given by the analytic expres- 
sion nini 



'LDA 



(k, W ) 



\bh z (ftV - E 2 re ) 



(ft V - E 2 r 



2E r 



Here, E rP = 



(6) 

is the recoil energy and fi = 
0.5hu; z (15\ 2 Na/a z ) 2 / 5 is the chemical potential. This 
analytic expression can also be used for 3D cigar-shaped 
trapped Bose systems. The above expression can be re- 
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where C = (15X 2 Na/a z ) 1 ^ 5 ~ 5.93 for the experiment Q 
and Cj — uj/u) z . In Fig. 2 we plot the dynamic structure 
factor S(k ) u>) vs the frequency oj by using the LDA as 
well as by using the Fourier transformation of the density- 
density correlation function given in Eq. ©. For finite- 
energy resolution we have replaced the delta function in 
Eq. © by the Lorentzian with a width of T = 0.1lo z . 



S (k, co) 
500- 



400 
300 
200 
100 



/ 



/ 




10 



co/co z 



FIG. 2: (Color online) Plots of the dynamic structure factor 
S(k,u>) vs normalized frequency u/uj z at T = 0. The dashed 
lines shows the dynamic structure factor based on the LDA. 

As one can see from Fig. 2, the dynamic structure factor 
has multiple peaks. This phenomena is due to the under- 
lying discrete axial spectrum. Fig. 2 shows that for the 
given k = 8, the most prominent peaks at co = Z.\Q2co z 
which corresponds to j = 5 and other small peaks cor- 
responds to the other quantum numbers. In order to 
observe the j ' = 5 quasiparticle modes, the wave num- 
ber of the Bragg pulse should be k — 8.0 which can be 
obtained by suitable choices of the photon wave number 
(k p ) and the angle (0) between the Bragg beams. Note 
that, k = 8.0 implies ka$ ~ 0.45 > k c ao. Therefore, 
the excited states will be separated out from the con- 
densate and the momentum measurements can be per- 
formed. The simplified LDA picture fails to produce the 
discrete structure of the axial modes, and it describes 
poorly the envelope of the spectrum of the dynamical 
structure factor. However, the single peak in Slda oc- 
curs at CjJ m ~ 4.75 for k = 8, but the locations of the 
highest peak in S(k,u>) is co/co z = 3.162. Therefore, the 
highest peak in S(k,uj) do not match with the peak in 
the SLDA(k,co). The LDA fails to describe the discrete 
structure of the modes because it assumes the system is 
locally uniform and it does consider partially the effect 
of the finite axial size. We have also carried out the same 
analysis for higher values of k. For k = 9, the condensate 
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responds resonantly at the frequency ui = 3.674w z which 
corresponds to the axial quantum number j = 6. In both 
the cases, the resonantly excited states are in the phonon 

regime since fc£ << 1, where £ = J h 2 /2m[i ~ 0.168a z is 
the healing length. 

The behavior of these multiple peaks in the dynamic 
structure factor can be resolved in two-photon Bragg 
spectroscopy, as shown by Steinhauer et al. In the 
two-photon Bragg spectroscopy, the dynamic structure 
factor can not be measured directly. Actually, the ob- 
servable in the Bragg scattering experiments is the mo- 
mentum transferred to the condensate which is related 
to the dynamic structure factor and reflects the behav- 
ior of the quasiparticle energy spectrum. The popula- 
tions in the quasiparticle states can be controlled by 
using the two-photon Bragg pulse. When the conden- 
sate is irradiated by an external moving optical potential 
V op — VB(t)cos(kz — tot), the excited states are popu- 
lated by the quasiparticle with energy Tilo and the mo- 
mentum hk, depending on the value of k and u) of the 
optical potential V op . Here, Vb is the intensity of the 
Bragg pulse. Suppose the system is subjected to a time- 
dependent Bragg pulse which is switched on at time t > 
and k is also along the z-direction. The momentum trans- 
fer to the Bose system from the optical potential can 
be calculated analytically either by using the Bogoliubov 
transformation 18] or by using the phase-density repre- 
sentations of the bosonic order parameter |l9l| and it is 
given by 



P z {t) = J^hk < a ](t)&j(t) >= 



j,k 



2h 



3 



x i ? J -[(w i -w),i]-i7 i [(w i +£<;),«], 



hkSj(k) 
(8) 



where ctj(t) is the time-evolution of the quasiparti- 
cle operator of energy Tiujj and Fj[(u>j ± oj),t] — 



n[(uj j ±u)t/2] 



For positive lo and a given k such that 



Sj(k) is maximum, the momentum transferred P z (t) is 
resonant at the frequencies lo — tOj. The width of the 
each peak goes like 2ir/t. In order to resolve the different 
peaks, the duration of the Bragg pulses should be at least 
of the order of the axial trapping period T z — 2ir/w z . 
Moreover, the duration of the Bragg pulse, is, (atom- 
light interaction time) is a main factor in the Bragg spec- 
troscopy and it must be of order ts > ni/hk 2 in order to 
populate the Bragg reflected beam significantly [i^J] ■ For 
large t and uj z « lo , P z (t) ~ S(k,ui) 0. In Fig. 3, we 
plot the net momentum transfer P z (t) vs the frequency 
lo for three different choices of the time duration of the 
Bragg pulse. Fig. 3 shows that the shape of the P z (t) 
strongly depends on the time duration of the Bragg pulse 
ts- When ts = 1-1T Z , the P z (t) is a smooth curve with a 
single peak at lo/lo z — 1.8, where T z = 2tt/lo z is the axial 
trapping period which is 40 msec for the experiment 
When ts = 1-9T Z , there is a little evidence of few small 



peaks start developing in the P z (t) . When ts = 3T Z , the 
multiple peaks in the P z (t) appears prominently. The 
location of the peaks in Fig. 3 for ts = 3T 2 are exactly 
same as in Fig. 2. It implies that P z (t) ~ S((k,ui) for a 
long duration of the Bragg pulse. Therefore, the multiple 
peaks in S(k, uj) are resolved in Fig. 3 only when the du- 
ration of the Bragg pulse is ts » T z . It should be noted 
that it is very difficult to use long duration of the Bragg 
pulse at the present situation. It can induce a dipole 
oscillation of the whole condensate in the trap |2lJ and 
the reflection of phonons at the boundaries which might 
causes a broadening of the response. We have checked 
that the P z (t) reflects the dynamic structure factor cal- 
culated from the local density approximation if is <T Z . 
We also found that Fig. 3 and the locations of the discrete 
peaks does not depend on the intensity Vb of the optical 
potential. We have also studied P z (t) for k = 9 and find 




FIG. 3: (Color online) Plots of the net momentum transferred 
P z {t) vs normalized frequency uo/u) z at T — when the di- 
mensionless wave-vector is k = 8 for various time duration of 
the Bragg pulse: ts = 1-1T Z (dashed), ts = 1.9T, (solid) and 
ts = 3.0T Z (dotted). Also, we have assumed Vb = 0.05?laj z . 

that there is a large peak at lj — 3.674cj z corresponds 
to the quantum number j = 6 if the time duration of 
the Bragg pulse is longer than the axial trapping period. 
This peak is also occurred in the dynamic structure factor 
as we have already discussed. 

One can do the same analysis for other values of k 
as long as only the low-energy axial excited states con- 
tribute significantly in the dynamic structure factor and 
these low-energy axial states should be less than the first 
radial excitation. It would be interesting to study in- 
termediate regimes of the energy and wave- vector where 
both the radial and axial discrete states are excited si- 
multaneously and it would provide a rich physics of the 
dynamic structure factor as well as the momentum trans- 
fer to the Bose system, due to the non-trivial structure of 
spectrum of the density fluctuations and its correspond- 
ing eigenfunctions. 

In conclusion, we have proposed an experiment in 
which the long wavelength discretized axial modes in a 
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cigar-shaped (but 3D) condensate can be imprinted by 
using the Bragg scattering experiments, similar to the 
recent experiment [j| where the short wavelength axial 
phonons with discrete radial modes have been observed. 
We have estimated the two-body scattering length a, also 
the values of the momentum (Kk) and energy (hu>) of the 
two-photon Bragg pulses which are needed in the Bragg 
scattering experiments in order to observe the low energy 
discrete axial modes. These discrete axial modes can be 
observed when the system is dilute enough and the time 
duration of the Bragg pulse is long i.e. ts >> T z . 
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